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Abstract 

We consider the realization of a quantum computer in a chain of nuclear spins cou- 
pled by an Ising interaction. Quantum algorithms can be performed with the help of 
appropriate radio-frequency pulses. In addition to the standard nearest-neighbor Ising 
coupling, we also allow for a second neighbor coupling. It is shown, how to apply the 
2'Kk method in this more general setting, where the additional coupling eventually allows 
to save a few pulses. We illustrate our results with two numerical simulations: the Shor 
prime factorization of the number 4 and the teleportation of a qubit along a chain of 
3 qubits. In both cases, the optimal Rabi frequency (to suppress non-resonant effects) 
depends primarily on the strength of the second neighbor interaction. 

PACS: 03.67.-a, 03.67.Lx, 03.65.Ta, 03.67.Dd, 03.67.Hk 

1 Introduction 

The Ising-spin chain has been proposed in (Lloyd 1993, Berman et al. 1994, Lloyd 1995, 
Berman et al. 2001) as a theoretical system which allows to implement a quantum computer. 
Typically, one would think of a chain of spin-1/2 nucleons embedded into a solid crystal, as 
a possible physical system, whose dynamics may be well described by the Ising Hamiltonian. 
This system must be subjected to a magnetic field, constant in time, with a sufficently strong 
variation along the spin chain. Additional RF-pulses (radio-frequency pulses) then allow the 
coherent control of the state of the system such that a quantum protocol can be realized 
(Berman et al. 2000, Berman et al. 2001, Berman et al. 2002). This is a special form of 
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NMR quantum computation as described in (Jones 2000). Ultracold atoms in optical lattices 
may provide an alternative physical realization of a Ising-Hamiltonian (Garcia-Ripoll and 
Cirac 2003). 

Typically, in these system we find two types of dipolar interactions: (i) Intrinsic dipole 
couplings between the nuclear spins. Those scale as distance and can be cancelled by 
the "magic angle" method (Slichter 1996). (ii) Mediated (mainly by the electrons) dipole 
couplings which have no clear distance dependence. 

Up to now, this model has been developed just theoretically and hopefully the technological 
and experimental part may start in a near future. However, because the Hamiltonian of this 
system is well known, many theoretical studies have been made (Berman et al. 1994, Berman 
et al. 2000, Berman et al. 2001, Berman et al. 2002^'^), Berman et al. 2002^^), Berman et 
al. 2002*^'^'', Lopez et al. 2003, Celardo et al. 2005) which are also important for the general 
understanding of quantum computation. In this model, first neighbor Ising interaction among 
the nuclear spins of paramagnetic particles of spin one half was considered. Thus, in this 
paper we want to consider also second neighbor Ising interaction among the nuclear spins. 
The transverse coupling will be neglected since one could expect that their coupling constants 
to be a least two or three orders of magnitude smaller than the longitudinal coupling constant 
(Ising) for our particular configuration. With a register of 4 qubits, we perform a numerical 
simulation of Shor's factorization algorithm (Shor 1994) of the number 4 and teleportation 
(Bennett et al. 1993) of an arbitrary qubit in a chain of three qubits, allowing for an interaction 
between second neighbors in the system. 



We consider an Ising spin chain with nearest and next-nearest neighbor interaction as a model 
for a quantum register. The spin chain is subject to a constant external magnetic field in 
^-direction, as well as to RF-pulses (with the magnetic field vector in the x-y plane). This 
chain is inside a strong magnetic field in the z-direction and may be subject to RF-pulses 
(with the magnetic field vector in the x-y plane). The constant magnetic field B[z), which 
must be extremely strong, also has a field gradient in the z-direction, which allows individual 
addressability of the qubits. During an RF-pulse, the whole external field may be written as 



where Bq, w and (p are the amplitude, the angular frequency and the phase of the RF-field. 
They are assumed to remain constant during a pulse, but are typically chosen differently for 
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The model 



B = {Bq cos{wt + (p), —Bq sin{wt + (p), B{z)) , 
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different pulses. Without any RF-field, the Hamiltonian reads: 

(n n— 1 n—2 \ 

J] I/;, 4^ + 2J 5^4^4^1 + 2/ 5^ 4^ 4^2 , (2) 
fc=l k=l k=l J 

where Wk is the Larmor frequency of spin k. We denote with |Ofc) the state where the nuclear 
spin k is parallel to the magnetic field and where it is anti-parallel. The RF-field induces 
the desired transitions between the Zeeman levels of the systems. 

The structure of the Hilbert space of the spin chain is particularly appropriate for quan- 
tum information studies, where the basis unit of information is a two-level quantum system 
("quantum bit", or qubit for short). Any such state \E' = Co|0) + Ci\l) can be represented 
with respect to some basis stats |0) and |1) by two complex numbers Cq and Ci such that 
|CoP + IC'iP = 1- The L-tensorial product of L-basic qubits form an L-register of L-qubits. 
In this space, we denote the resulting product basis by \a) = \iL-i, ■ ■ ■ ,io) with ij = 0, 1 for 
j = 0, . . . , L — 1. A pure wave function can be expanded in this basis by \1/ = ^ Ca\a), where 
J2a l^aP — For notational convenience, we require that a = X]j'=o^ h '^^ ■ 

The Hamiltonian Hq in Eq. (2) is diagonal in the computational basis (the product basis 
defined above): 

(-l)°fc 

Hq \an-i ■ ■ ■ «iao) = Ea \an-i ■ ■ ■ a2ao) IlWk) = — - — |afc) 

, /n— 1 n—2 n—3 \ 

K = - I $^(-1)"'= + J^(-l)"'=+"'=+^ + J' ^(_i)«.+"'=+2 . (3) 

\fc=0 k=Q k=Q J 

The index a without subscript denotes the positive integer represented by the string a^-i • • • aiao 
in the binary number system, e.g.: Hq |101) = 1 101). Choosing the Larmor frequencies 
Wk such that Wk : Wk/wk-i = 2, leads to a spectrum which (ignoring the spin-spin coupling 
terms) has equidistant levels: — E^-i = constant. 

The RF-pulses are essential for any implementation of a quantum algorithm. During such 
a pulse, the full Hamiltonian may be written in the form 

H = Ho + W{t) W{t) = Yl (e'^"*^''^ + e-'^"'*^^) 4") , (4) 

k=l 

where the frequency w, the phase-offset ip and the Rabi frequency Q are free parameters 
(f2/7 is the amplitude of the RF field, where 7 is the nuclear spin gyromagnetic ratio). In 
the computational basis |«n_i . . .aiao) the raising and lowering operators for spin k can be 
written as: = |Ofc)(lfc| and = |lfc)(Ofc|. We will consider sequences of pulses, where each 
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pulse may vary in time. During the pulses and Vt are fixed as well as the RF frequency w 
which is assumed to be on resonance with some allowed transition. We assume that there are 
no degeneracies. Hence, 



E, 




Wk + J[ (-1)"^+! + ]+J'[ (-1)°'=+^ + (-1)"^-^ ] , (5) 



w = 



where it is understood that (—1)"' = if / < or / > n. Thus, in the center of the 
spin chain (away from the borders) the J-coupling may lead to frequency shifts of Aw = 
or ±2 J {Aw' = or ±2 J'), whereas at the borders the frequency shifts are Aw = ±J 
(Aw' = ±J'). The following notation is based on these observations. The unitary evolution 
(in the interaction picture) during a resonant RF pulse is denoted by 



with /i, z/ G [—2,-1,0,1,2]. Note that in order that R'^''^{QT,ip) is a resonant pulse, the 
indices fi, u, k must fulfill certain relations as discussed above {i.e. not all combinations lead 
to resonant transitions). The wave function dynamics during such RF-pulses is computed 
numerically as sketched in App. A. 

In the numerical simulations to follow (Sees. 4 and 5), we have chosen the following 
parameters in units of 27rxMhz, 



These parameters give rise to a simple Zeeman spectrum with equidistant levels. For the 
possibly more realistic design described in (Berman et al. 2002*^^^), we would expect similar 
results. In our simulations, we will take advantage of the second neighbor interaction in order 
to reduce the number of pulses for the realization of the quantum algorithms (Lopez and Lara 
2006). Those are Shor's factorization in Sec. 4 and teleportation in Sec. 5. 

3 Second neighbor interaction and the 27rfc-method 

For the considerations in the previous section, only resonant transitions have been taken into 
account. That allowed trivially to describe the dynamics under the RF-pulses analytically. 
However, it is possible to go beyond this simple picture. To that end we distinguish near 
resonant transitions where the frequency differs from the resonant frequency by values of the 
order of J or J' and far resonant transitions, where the difference is of the order of the Larmor 
frequencies. Then it can be shown that the Hamiltonian decomposes into independent 2x2 
matrix blocks as long as far resonant transitions are neglected. This still allows to describe 




W = Wk + fiJ + uJ' , 



(6) 



ujo = 100 , uji = 200 , UJ2 = 400 , us = 800, J = 10 , J' = 0.4 , = 0.1 . 



(7) 
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the dynamics under the RF-pulses analytically, and in particular it allows to control and 
suppress the non (near) resonant transitions. The 27r/c-method (Berman et al. 2002*^") ) and 
its generalization (Berman et al. 2002*^^^^) resulted from such considerations. 

Assume we perform the pulse Ri'~^{flT,ip) in the 4-qubit quantum register. This pulse 
induces the resonant transitions 



at the frequency w 
frequencies w + A 



|0001) ^ lOOll) lOlOO) ^ lOllO) (8) 
Wi — J'. However, it also induces the near-resonant transitions at 



|1001) ^ |1011) 



|1100) 


!> 


|1110) 


A = 


2J' 


loooo) 




lOOlO) 


A = 


2J 


loioi) 


<— > 


lOlll) 


A = 


-2J 


|1000) 




|1010) 


A = 


2J + 2J' 


|1101) 




lllll) 


A = 


-2 J + 2 J' 



(9) 

In the near resonant approximation, we may write down an evolution equation for each 
transition separately. To this end let a denote the first state of the transition pair (where 
«! = 0) and /3 the second (where /3i = 1, V/c 7^ 1 : /3k = Ofc)- Then we find for the wave 
function coefficients in the interaction picture (see App. A): 



The solution for the evolution operator is 

= U{t) 



(10) 





cos%^ - i| sin%^ 
g e"'"^ sin %^ 



g e*^ sin %^ 
cos + sin ^ 



(11) 



where fie = Vfi^ + A^. Hence, e.g. for a 7r-pulse t = r = 7r/f2 at the end of the pulse, and we 



may tune Q such that any near resonant transition is switched off. This only requires that 

TT Q, 



2n 



A2 



(12) 



The 27rfc-method as presented here has two shortcomings. If there are several near-resonant 
transitions with different frequency shifts A it will be impossible, in general, to eliminate all 
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of them. Second, any near resonant transition also implies a phase rotation, which cannot 
be corrected by the method. In principle some of those shortcomings can be overcome at the 
expense of more complex pulse sequences, as shown in (Berman et al. 2002c). 

For the numerical simulations in Sees. 4 and 5 we define the optimal Rabi frequency to 
eliminate a near-resonant transition with shift A during a vr-pulse as 

fij) = . (13) 

Note that for a 7r/2-pulse the corresponding optimal Rabi frequency is given by the same 
equation, but with k replaced by 2k. 



4 Simulations of Shor's factorization algorithm of num- 
ber four 

An experimental realization of Shor's factorization algorithm has been demonstrated recently 
in (Vandersypen et al. 2001) using nuclear magnetic resonance. Following Shor's approach 
(Shor 1994) for factorizing an integer number A^, one selects a L + M-register of the form 
|x; y), where is the input register of length L, and \y) is the evaluation register of length 
M. The y-register is used to store the values of the periodic function y{x) = g^(mod A^), 
where the integer q is chosen co-prime to A^ [i.e. the greatest common divisor gcd(g, A^) = 1). 
The algorithm is divided into three parts (Nielsen and Chuanf 2000, chapter 5). First, one 
creates the uniform superposition of all basis states in the x-register. Second, one chooses 
some co-prime number q and computes the function y{x) = g^(mod A^) in the y- register. 
Third, one applies the inverse discrete Fourier transform to the x-register. After these steps, 
one measures the state in the x-register, which provides the information on the factors of A^. 
For more details, we refer the interested reader to (Nielsen and Chuang 2000). 

For factorizing the number A^ = 4, two qubits in each register are sufficient {L = M = 2). 
The only co-prime number of A^ is g = 3, and thus the period of the function y{x) = 3^(mod 4) 
gives the factors of A^. In the present case this period is T = 2. Applying the procedure above 
to the initial state \I^o = 1 00; 00), one obtains after each of the three steps: 

^(ideal) ^ 1 ^ ^ ^ ^ ^ 

^(idCal) ^ 1 ^ ^ ^ I^Q. ^ 1^^. ^ ^^^^ 

^(ideal) ^ 1 ^ ^ ^ I^Q. ^ ^^^^ ^ ^ ^^^^ 
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The measurement on the x-register give us the states |00) or |10) (x = or x = 2), which 
imphes that the period of the function y{x) is T = 2, as expected. 

The following computation is performed in the interaction picture. Using the parameters 
from Eq. (7) and starting with the ground state of the system, \l/o = |0000), we create 
with the help of three 7r/2-pulses: 

v&i = R^'^\n/2, n/2) R}^\tx/2, n/2) RY{n/2, tx/2) M/q . (17) 

The evaluation of the function y{x) = 3^' (mod 4) in the y-register is carried out with four 
TT-pulses: 

^2 = Rf'^'in, 7r/2) Rf'\7r, n/2) Rl'~\7r, tx/2) R]^\tx, tx/2) . (18) 
Finally, the discrete Fourier transformation in the x-register is obtained through five vr-pulses: 

^3 = /22~''"'(^,^/2) i?o''"'(7r,-7r/2) K^'-\tx,tx /2) R''/{tx, -tx/2) R,'^^\tx, -tx/2) ^2 ■ 

(19) 

The final measurement consists in tracing out the y-register, and measuring the probabilities 
of finding the x-register in any of the possible basis states. In the present case, we would get 

(00|f?,|00) = (10|f?,.|10) = ^ (OII&IOI) = (ll|f?.|ll) = , (20) 

where = try|\E'3) (\E'3|. This yields the expected period T = 2. 

In what follows, we simulate the pulse sequence for Shor's algorithm with our model 
Hamiltonian of Eq. (4). The wave-packet evolution is carried out in the interaction picture 
using the parameters given in Eq. (7). In Fig. 1 we show the difference between the final state 
\E'3, obtained from the evolution with the full Hamiltonian H{t), and the ideal final state 
^^ideai)^ The latter is obtained from an evolution where only the resonant transitions have 
been taken into account. As deviation, we plot the real (red impulses) and imaginary (blue 
impulses) parts of the difference between the respective expansion coefficients Da and Da^'^'^'^^ 
as a function of a, which denotes the basis states, as explained below Eq. (3). We can see 
that by far the largest error consist in an accumulation of residual phases in the expansion 
coefficients of those basis stats which are expected to be populated in \[f 3 , i.e. |0001), 
|0011), |1001), and |1011). This can be deduced from the large imaginary parts errors, visible 
in Fig. 1. However, the probabilities of the states remain within their right value. 

In order to quantify the deviation of the real evolution with H{t) from the ideal one, where 
only resonant transitions are taken into account, we introduce a figure of merit, the fidelity 
(Peres 1984, Schumacher 1995, Gorin et al. 2006) 

^ _ |^^(idcal)|^^|2 ^ (21) 
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Figure 1: The difference in the expansion coefficients between the final state \l/3 obtained 
from the evolution with the full Hamiltonian and the ideal final state vl/^''^°^'''. Red impulses 
(towards the left) show the real part, blue impulses (towards the right) the imaginary part. 
The difference in the expansion coefficients is plotted versus a, which labels the basis states. 



8 



0.65 ^ ^ ^ ^ ^ ^ ^ ^ 1 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

Q 

Figure 2: The fidelity, Eq. (21), for tlie final state resulting from the factorization algorithm 
as a function of the Rabi frequency Q. The blue arrows show the optimal Rabi frequencies 
Jl^'* according to Eq. (13), for A = 2J' and k = 1,2, 3, and 4 (from right to left). 

which gives the probability to find the system in the state even though it really is in the 

state As we have seen in Sec. 3 that one can influence the accuracy of the implementation 
of the quantum protocol by appropriately choosing the Rabi frequency (the strength of the 
RF field) we will first investigate the dependence of the fidelity -Fg = |(\E'3'^'''^^''|\i/3)p on 
the Rabi frequency. The result is shown in Fig. 2, where is plotted versus Q in the 
interval (0.08, 0.48). The blue arrows show the positions of the optimal Rabi frequencies 
according to Eq. (13) for A = 2 J' and k = 1,2,3, and 4 (from right to left). The fidelity 
clearly behaves as expected. We find maxima for those values, where the 27r /c-relation is 
fulfilled. The resulting curve is rather smooth at small Rabi frequencies, while additional fast 
oscillations appear near Q ^ 0.3. The additional structures are due to the transitions where 
the detuning is of the order of 2 J. Since J ^ J' these transitions have nearly no effect for 
small Rabi frequencies, but there effect becomes more and more significant as VL increases. In 
practice, one would have to find a compromise between a sufficiently fast evolution (for that 
one needs large Rabi frequencies) and a sufficiently high fidelity, which can be achieved more 
easily at small Rabi frequencies. 

In order to study the mechanism behind the 2ti A;-method in more detail, we record the 
(loss of) fidelity during the execution of Shor's quantum algorithm as a function of time. The 
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0.93 I \ ^ ^ ^ ^ ^ ^ 1 

-50 50 100 150 200 250 300 350 

t , t - 50 

Figure 3: The fidelity, Eq. (21), during the quantum factorization algorithm as a function of 
time t. We show two cases: ^2 = = 0.1 for A = 2J' (red and blue lines), and Q = 0.116 
(light blue and pink lines), where the fidelity of the final state has a minimum (see Fig. 2). 
For each case, we used alternating colors to distinguish the individual pulses in the pulse 
sequence Eqs. (17-19). For better visibility, the curve for Q = 0.116 is shifted to the left by 
50 time units. 
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result is shown for two cases, n = n^^> = 0.1 for A = 2J', and Q = 0.116 in Fig. 3. In the first 
case, the Rabi frequency is a optimal one, according to Eq. (13), whereas in the second case, 
the Rabi frequency is such that Fg, as plotted in Fig. 2 has a minimum. In both cases, we see 
the expected Rabi oscillations (in F{t)) during the majority of the vr-pulses. In the first case, 
the oscillations are completed during the vr-pulses, such that F(t) at the end of a vr-pulse is 
almost returns to the value it had at the beginning of the pulse. By contrast, in the second 
case, the vr-pulses stop when F(t) is almost exactly in a minimum of the near-resonant Rabi 
oscillations. In the graph of the second case (light blue and pink lines) we can see that for 
some pulses {e.g. the 8'th and the ll'th) there are apparently no oscillations, only a striking 
broadening of the curve for F{t). During those pulses, the near-resonant transitions with 
A = 2 J' involve states which almost unoccupied, such that the transitions cannot really take 
place. In those instances we have a chance to see the effects of transitions with A being of 
the order of 2J. The resulting oscillations in F{t) are too fast to be resolved, such that the 
only visible effect is a broadening of the curve. This also demonstrates that the transitions 
with detuning of the order of 2 J can be neglected at small Rabi frequencies. That explains 
the smooth behavior of Ft[{Q), depicted in Fig. 2. 



5 Quantum teleportation on 3-qubits 

The basic idea of quantum teleportation (Bennett et al. 1993) is that Alice (left end qubit 
in our chain of three qubits) and Bob (the other end qubit) share two qubits which are in 
a maximally entangled (Bell) state. Here, we shall follow the prescription of (Nielsen and 
Chuang 2000). 

$e = 4H (|0a00b) + |1a01b)) . (22) 
V2 

We adjoin to Alice the arbitrary state 

$. = ColO)+Cf|l) , (23) 

to be "teleported" to Bob. This results in the quantum state ($i = $x ® $c) of the whole 
system 

V 2 

Then, we apply a controlled-Not (CN) operation between the added qubit and that of Alice 

CN32 |23, 22, io) = Ks, ^2 © «2, «o) , (25) 



11 



where ^2 © ^3 = (^2 + i3)mod 2. This results in the state 



(ideal) 



[CI 10000) + Q 10101) + |1100) + Ci" 1 1001)) 

v2 



(26) 



In (Nielsen and Chuang 2000) a final Hadamard gate is applied to the added qubit. It allows 
Alice to measure the two qubits in the computational basis instead of the Bell basis; see 
(Bennett et al. 1993). In our implementation, we replace the Hadamard gate with a different 
7r/2-qubit rotation A, which however achieves the same goal: 



A, 



f|0000) _ 1 Jioooo) - |1000) 
||iooo) ~ 71 [loooo) + |1000) 



(27) 



On the Bloch-sphere, this corresponds to a rotation about the angle — 7r/2 around the j/-axis. 

>adcal) ^ ^(ideal) 



After some rearrangements, the final state $3 



reads: 



(ideal) 



- (100) ® |0) ( CI |0) + CI |1) ) + |01) ® |0) ® ( CI 1 1) + CI |0) ) 
+ |10) ® |0) ®{Cl\\)- CI |0) ) + 111) ® |0) ® ( CI |0) - |1) )} 



(28) 



When Alice measures both of her qubits, there are four possible cases: |00), |01), |10), and 
1 11). For each case. Bob will get the original state <I>x, provided he applies the proper operation 
to his qubits. The operation he has to choose depends on the outcome of Alice's measurement 
as follows: 



Alice's result 


Bob's operation 


|00) 


id 


|01) 


N 


|10) 




111) 





where 



id 

N 

(Jr. 



|0) 
|0) 
|0) 



|o), |i) 

ll), ll) 

|o), |i) 



ll) 

|0) 



(29) 



|1) 



Note that neither Alice nor Bob need to know the state $x in order to perform the transfer 
to Bob's qubit. In fact if any of the two "knew" something about the state, the state transfer 
would not work. 

In order to implement the quantum teleportation scheme on our Ising spin chain quantum 
computer, we start with a spin chain of 4 spins (0, ... 3), where spin 3 is in the unknown state 

$0 = d lOOOO) + CI 1 1000) . (30) 
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Figure 4: The difference in the expansion coefficients between the final state $3 obtained 



from the evolution with the full Hamiltonian and the ideal final state $ 



(ideal) 



Red impulses 



(towards the left) show the real part, blue impulses (towards the right) the imaginary part. 
The difference in the expansion coefficients is plotted versus a, which labels the basis states. 



The first step consists in generating an entangled (Bell) state between qubit 2 (Alice) and 
qubit (Bob). This is obtained with the following three pulses: 

$1 = Rl^~\n, -n/2) R'/{n/2, -n/2) Rl'\n/2, -7r/2) $0 • (31) 

Then the controUed-not operation CN32 is applied through 

$2 = i?°-^(7r, -n/2) R^/{7i, n/2) $1 . (32) 

The final state, given in Eq. (28), is obtained after applying the A-^ to qubit 3 via the following 
two pulses: 

$3 = Rs''\n/2, -7r/2) R'/iiT/2, -7r/2) $2 . (33) 

As an example, we choose Cg = 1/3 and Cf = V^/S for the coefficients of the unknown state. 

Again, we simulate the pulse sequence with our model Hamiltonian of Eq. (4) using the 
parameters of Eq. (7), just as in the previous section. Fig. 4 shows the difference between the 
final state $3, obtained from the simulation and the ideal final state $3'^°'^''', given in Eq. (28). 
The figure is produced exactly in the same way, as Fig. 1 for the case of Shor factorization. 
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0.65 ^ ^ ^ ^ ^ ^ ^ ^ 1 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

Q 

Figure 5: The fidelity, Eq. (21), for tlie final state resulting from the teleportation algorithm 
as a function of the Rabi frequency Q. The blue arrows show the optimal Rabi frequencies 
Jl^'* according to Eq. (13), for A = 2J' and k = 1,2, 3, and 4 (from right to left). 

Here, we find somewhat smaller deviations than in Fig. 1, which may be simply explained by 
the number of pulses for the teleportation protocol being smaller. 

In Fig. 5 we show the fidelity = "''^^|<l>3)|2 function of the Rabi frequency 
Q. Again, Fig. 5 is the precise analog of Fig. 2, where the same quantity is plotted for the 
Shor factorization protocol. We again find the expected maxima at those points, where the 
27r fc-relation, Eq. (13), is fulfilled. At least for sufficiently small frequencies, there are no 
noticeable effects of the nearest neighbor interaction, which implies much larger detuning 
A ~ 2 J. The present protocol contains vr-pulses, as well as 7r/2-pulses. However, note that 
for 7r/2-pulses, the 27r fc-relation is fulfilled for even k, only - see discussion below Eq. (13). 
On Fig. 5, the corresponding frequencies are marked by the second and the forth arrow, 
respectively (counting from the right). The data does not show any modulation with respect 
to k being odd or even. This is a clear signature that the fidelity of the present protocol is 
largely determined by the vr-pulses alone. 

Fig. 6 shows the (loss of) fidelity during the execution of the teleportation protocol as a 

(2) 

function of time. In the present case, we choose n = n'^^' = 0.1 for A = 2J'. Thus, the Rabi 
frequency fulfills the 27r A;-relation, Eq. (13), such that the final fidelity is maximal, locally (see 
Fig. 5). In Fig. 6 we can distinguish the 7i/2 pulses (two on each side) surrounding the three 
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Figure 6: The fidelity, Eq. (21), during the quantum factorization algorithm as a function of 
time t for Q = =0.1, where A = 2 J' (red and blue lines). We used alternating colors to 
distinguish the individual pulses in the pulse sequence Eqs. (31-33). 
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TT-pulses in the middle of the pulse sequence. Oscillations are visible only during the central 
vr-pulses, not during any of the 7r/2-pulses. Note, however that during the last 7r/2-pulse we 
find very fast oscillations with a relatively small amplitude, which must come from detuning 
of the order of A ~ 2J. The behavior of F[t) clearly demonstrates that the loss of fidelity 
occurs almost exclusively during the execution of the 7r-pulses. 

6 Conclusion 

For a quantum computer realized with a one-dimensional chain of nuclear spins (one half), 
we have studied the effect of an additional second neighbor Ising interaction. This allows 
greater flexibility in the choice of appropriate pulse sequences for the implementation of a 
given quantum algorithm. We have found that after adapting the 27r/c-method to this new 
situation, the desired gate operations can be realized with high fidelity. We have illustrated 
our results with two case studies: Shor's quantum algorithm for the factorization of an integer 
number, and another algorithm which allows the teleportation of a qubit across the spin chain. 
In both cases, we have studied the fidelity of the whole algorithm as a function of the Rabi 
frequency, as well as the decay of fidelity during the algorithm as a function of time. In spite 
of the fact that both algorithms contain vr- and 7r/2-pulses, we found that only the 7r-pulses 
were responsible for the decay of fidelity. We found also that as long as the nearest-neighbor 
interaction J is much larger than that of the next-nearest neighbor interaction J', the (loss 
of) fidelity is dominated by the latter. In other words, if the second neighbor interaction in 
a chain of nuclear spins is not entirely negligible, it will typically dominate the non-resonant 
effects. 
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A Numerical solution of the Schrodinger equation 

We consider the time evolution of the quantum state \l/(t) of the Ising spin chain during a RF 
pulse as defined in Eq. (4). The Schrodinger equation reads: 

ihdt ^{t) = Hit) ^it) H{t) =Ho- — Y, (e'^"*^''^ It + e-^("*+^) J,-) , (34) 

fc=i 
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where \I/(0) is the quantum state at the beginning of the pulse, and \I/(r) is the state at the 
end (r is the duration of the pulse). The unperturbed Hamiltonian Hq is given in Eq. (2). 
Expanding in the basis of the time-independent part Hq, we arrive at the matrix form 
of the Schrodinger equation: 



In order to obtain the dynamics in the interaction picture, we set 

2 1 e-'l'^>'+«'lD„i„...„(() : a,. = 1 



fc=0 I 

where 



^k = w-Wk-J[{-lT''+^ + {-lT'-'] - J' + , (37) 

as obtained from Eq. (5). Here, it is understood that (— 1)°' = if / < or / > n. In order 
to evolve the wave packet numerically, it is convenient to rewrite the system of differential 
equations in terms of real variables. To that end, we introduce the real variables Xct{t) and 
Yait) such that = + 'lY^. Then we find: 

n— 1 ( 

d X {t) -—S^ { ^^^'^^'^'^ + Xa\a,,=l{t) - COs(Afct + (p) Ya\ak=i{t) : Ofc = 

* " 2 ^ |sin(Afct + V?) X„|c,^=o(^) - cos(Afet + v?) ya|afe=o(^) : "fc = 1 

d^Y {t) = -y^ |cos(Afet + ^) X„|«,=i(t) - sin(Afct + if) F„|«,=i(t) : afc = ^^^^ 
2 fcl^ + V') -^a|afc=o(^) + sin(Afct + V9) Fa|Q^=o(^) : «fc = 1 
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